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Abstract 

Teleoperation  provides  a  human-machine 
interface  that  allows  an  operator  to  control 
machines  at  a  large  distance.  It  has  recently 
been  the  subject  of  a  significant  amount  of 
research  in  the  robotics  community. 
Teleoperation  can  also  provide  an  instructor 
with  an  excellent  motivating  example  for  a 
control  systems  course  where  it  can  serve  as  a 
platform  for  introducing  several  key  controls 
issues  such  as  noise  compensation  and  time 
delay.  Computer  packages  such  as  Matlab  and 
SIMULINK  are  readily  available  in  most 
engineering  schools  for  simulating 
teleoperation  systems.  In  this  article, 
SIMULINK  is  used  to  present  a  case  study  of 
a  simple  teleoperation  system.  Initially,  a  PI 
controller  is  introduced  to  control  a  simple 
bilateral  teleoperator  model  without  noise  or 
time  delay.  It  is  then  shown  that  the  inclusion 
of  time  delay  can  significantly  affect  the 
system  performance  and  can  even  cause  the 
system  to  become  unstable.  A  relatively  new 
technique  called  the  wave  variable  method  is 
then  used  to  stabilize  the  system  for  any  fixed 
time  delay.  Lastly,  recommendations  for 
further  educational  projects  involving 
teleoperation  are  given. 

I.  Introduction 

Since  the  introduction  of  the  first  modem 
master/slave  manipulator  in  the  late  1940’s, 
teleoperation  systems  have  been  used  for  a 
number  of  different  tasks,  e.g.,  handling  toxic  or 
harmful  materials,  operating  in  remote  areas 


Figure  1:  A  bilateral  teleoperation  system. 


such  as  undersea  or  space,  and  performing  tasks 
that  require  extreme  precision,  and  will  continue 
to  play  an  increasingly  important  role  for  such 
applications  in  the  future  [5].  Recently, 
teleoperation  has  gained  greater  public  awareness 
through  the  Mars  Rover  project.  Future 
applications  include  robot  surgery  in  remote 
communities,  war  zones,  and  disaster  areas  and 
teleoperation  over  the  Internet.  There  are  also 
numerous  educational  applications  [12,  16]. 

During  teleoperation,  the  operator  typically 
uses  a  visual  display  and  a  “master”  manipulator 
(e.g.,  a  joystick)  to  manually  control  a  remote 
“slave”  device  such  as  a  robotic  arm.  The 
teleoperator  system  consists  of  a  human  operator, 
the  master,  the  communication  channel,  the 
slave,  and  the  environment.  To  improve  the  task 
performance,  the  contact  force  of  the  slave  with 
the  environment  can  be  reflected  back  to  the 
operator.  Adding  force  feedback  to  a 
teleoperation  system  provides  the  operator  with 
increased  awareness  and  can  considerably 
improve  his  or  her  ability  to  perform  complex 
tasks,  particularly  when  visual  information  is 
limited.  In  this  case,  the  system  is  called  a 
bilateral  teleoperation  system. 

Figure  1  is  an  illustration  of  a  haptic  bilateral 
teleoperation  system.  The  haptic,  or  force 
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feedback,  interface  provides  a  kinesthetic  link 
between  the  operator  and  a  virtual  environment 
giving  the  user  the  ability  to  actually  feel  the 
environment  surrounding  the  slave  manipulator. 
The  goal  is  for  the  slave  device  to  track  the 
master’s  behavior  with  respect  to  force,  position, 
and  velocity.  Under  ideal  situations,  this  goal  is 
easy  to  accomplish.  However,  when 
teleoperation  is  performed  over  a  large  distance, 
or  when  using  a  slow  communication  channel, 
time  delay  appears  in  the  information 
transmission  between  the  local  and  remote 
environments.  This  can  cause  a  serious 
degradation  in  system  performance.  Time  delay 
can  also  be  quite  irritating  to  the  user, 
particularly  when  the  user  is  forced  to  adopt  a 
“move-and-wait”  strategy,  whereby  a  small 
movement  is  made  and  the  operator  waits  to 
observe  the  results  of  the  movement  before 
committing  to  further  action.  Even  more  serious, 
time  delay  can  cause  a  force  feedback  system  to 
become  unstable.  System  noise  can  also  be 
problematic. 

These  issues  are  important  subjects  of  ongoing 
research  in  the  area  of  bilateral  teleoperator 
control.  Recent  results  in  the  area  are  promising, 
and  some  of  the  techniques  are  sufficiently 
developed  now  that  they  are  suitable  for 
inclusion  in  undergraduate  and  first  year 
graduate  courses  on  control  systems.  The  goal  of 
this  paper  is  to  describe  one  of  the  most 
promising  bilateral  teleoperation  control 
techniques  that  is  particularly  suitable  for 
educational  purposes.  While  a  significant  amount 
of  effort  and  expense  would  be  required  to 
construct  and  use  actual  teleoperation  hardware, 
the  material  proposed  in  this  article  only  requires 
software  models  that  can  be  easily  implemented 
using  readily  available  computer  packages  such 
as  Matlab  with  SlMULINK  [4,6,7,11], 
Throughout  the  paper  we  use  Matlab  and 
SIMULINK,  and  we  assume  that  the  reader  is 
familiar  with  these  software  packages;  however, 
other  similar  control  system  simulation  software 
can  be  used. 


The  next  section  introduces  a  simple  model  of 
a  bilateral  teleoperator  operating  under  ideal 
conditions.  In  this  case,  a  simple  PI  controller 
will  suffice.  After  simulating  the  ideal  case,  we 
introduce  the  student  to  some  “real  world 
problems”  such  as  noise  and  time  delay.  The 
problem  of  compensating  for  time  delay  in  the 
system  is  addressed  in  Section  III  by 
implementing  a  passivity-based  approach  called 
the  wave  variable  method.  Lastly,  conclusions 
and  recommendations  for  additional  exercises  to 
illustrate  teleoperator  control  issues  appear  in 
Section  IV. 

II.  Building  a  Matlab  and  SIMULINK  Model 

Initially  we  will  use  the  following  linear  model 
for  the  equations  of  motion  for  the  master  and 
slave  manipulators,  respectively, 

rc  =  J„0  +  Bs6 

where  xm,  xc  are  the  input  torques,  Jm,  Js  are  the 
inertia  terms,  and  Bm,  Bs  are  the  damping  terms 
of  the  system  being  simulated.  For  the 
simulations  shown  in  this  paper  the  following 
values  were  used:  Jm-Js=4  Kg-m2  and  Bm=Bs=2 
Kgm2/s.  To  determine  the  force  necessary  to 
move  the  slave,  the  following  simple  PI- 
controller  is  used: 

rc  =  Kpe(t)  +  K,e,  (t)  (2) 

where  e(t)  is  the  difference  between  the  desired 
slave  velocity  exJ(t)  and  the  actual  slave  velocity 
es(t),  ej(t)  is  the  integral  of  the  error,  and  Kp  and 

Kj  are  the  proportional  and  integral  gains  of  the 
controller,  respectively. 

The  SIMULINK  model  of  the  teleoperation 
system  requires  only  a  few  simple  blocks  as 
shown  in  Figure  2.  For  the  master  and  slave,  the 
transfer  function  block  is  used.  The  Pi-controller 
for  the  system  consists  of  several  simple  blocks. 
A  gain  block  for  the  proportional  gain  Kp  and 
another  gain  block  in  series  with  an  integrator 
block  for  the  integral  gain  Kj  are  added  together 
to  form  the  controller.  For  the  input  to  the 
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Figure  2:  SIMULINK  block  diagram  of  a  bilateral  teleoperation  system.  The  blocks  inside  the  dashed 
lines  represent  additive  noise  while  the  blocks  inside  the  dotted  lines  represent  time  delay. 
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Figure  3:  Torque,  position,  and  velocity  data  for 
the  simulated  model  of  a  bilateral  teleoperation 
system. 


Time(sec) 


Figure  4:  Torque,  position,  and  velocity  data  for 
the  simulated  model  of  a  bilateral  teleoperation 
system  with  a  time  delay  given  by  2T=500msec. 


system,  any  source  block  may  be  used;  for  the 
simulations  shown  a  step  function  block  was 
used  and  it  was  assumed  that  the  slave  was  in 
free  space  so  the  environment  input  was  zero. 
Since  the  outputs  of  the  master  and  slave  transfer 


function  blocks  are  velocities,  integrator  blocks 
must  be  used  to  obtain  position  information.  In 
order  to  have  access  to  the  output  data,  six  “to 
workspace”  blocks  were  used.  Each  block  was 
given  a  specific  variable  name  that  was  saved  in 
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the  Matlab  workspace.  Once  the  SIMULINK 
simulation  is  executed  and  the  output  data  is 
stored,  another  Matlab  routine  plots  the  data. 
This  file  reads  in  all  of  the  data  and  conveniently 
displays  the  variables  of  interest;  an  example  of 
this  is  shown  in  Figure  3.  In  this  case  we  have 
chosen  to  display  the  torque,  position,  and 
velocity  data  for  both  the  master  and  slave  in  a 
side-by-side  comparison.  Also,  the  data  can  be 
used  to  determine  other  variables  such  as  torque, 
position,  and  velocity  difference,  simply  by 
subtracting  the  slave  data  from  the  master  data. 

As  mentioned  earlier  once  a  problem  like  noise 
or  time  delay  is  added  to  a  system,  the 
performance  degrades  substantially.  In  fact,  if  the 
time  delay  is  sufficiently  large  the  system  will 
become  unstable.  Noise  can  be  added  by 
including  the  blocks  shown  inside  the  dashed 
lines  in  Figure  2.  In  this  article  we  will  not 
discuss  issues  concerning  the  control  of 
teleoperation  systems  in  the  presence  of  noise; 
however,  there  are  a  number  of  recent  papers  on 
this  important  problem  available  to  the  interested 
reader,  e.g.,  [1,13,17],  Instead,  this  article  will 
focus  on  time  delay.  Time  delay  can  be 
implemented  by  including  the  blocks  shown 
inside  the  dotted  lines  in  Figure  2.  For  all  of  the 
simulations  shown  in  this  article  we  assume  that 
Tr=Tl=T  where  Tr  is  the  delay  from  master  to 
slave  and  TL  is  the  delay  from  slave  to  master.  In 
this  case,  we  refer  to  the  total  round  trip  delay  as 
2T.  The  same  system  in  Figure  3  was  simulated 
with  2T=500msec,  and  the  results  are  shown  in 
Figure  4.  As  can  be  seen,  the  system  quickly 
became  unstable.  The  next  section  describes  a 
solution  to  the  problem  of  time  delay  by 
implementing  the  wave  variable  method. 

III.  Compensating  for  Instability  Due  to  Time 
Delay 

One  of  the  major  problems  that  can  be  found  in 
bilateral  teleoperation  systems  is  time  delay. 
Over  relatively  short  distances,  the  time  delay 
may  not  be  noticeable;  however,  when  the  master 
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Figure  5:  Pole-zero  plot  of  a  teleoperation  system 
with  a  time  delay  ranging  from  2T=0sec  to 
2T=10sec.  With  no  time  delay,  the  system  is 
stable,  but  as  the  time  delay  increases,  the  system 
becomes  unstable. 

and  slave  are  located  at  a  far  distance  from  each 
other,  time  delay  is  no  longer  negligible.  In  the 
case  of  space  teleoperation  the  roundtrip  delay 
can  approach  6sec,  and  surface  to  deep  ocean 
transmission  can  take  anywhere  from  seconds  to 
minutes  depending  on  depth  and  mode  of 
communication  [14].  Until  relatively  recently  [3], 
it  was  thought  that  time  delay  precluded  the  use 
of  force  feedback  in  such  cases. 

Time  delay  in  a  control  system  introduces  a 
phase  lag,  which  in  turn  degrades  the  system 
performance  and  can  cause  instability.  Figure  5 
shows  the  effect  that  time  delay  has  on  the 
location  of  the  poles  of  the  transfer  function  of 
the  original  system  described  in  Section  II.  The 
figure  was  generated  by  including  a  Pade 
approximation  of  the  time  delay  in  the  transfer 
function  of  the  overall  teleoperation  system  and 
using  a  ‘for’  loop  to  gradually  increase  the 
amount  of  time  delay  in  the  system  from 
2T=0sec  to  2T=10sec.  As  the  time  delay  is 
increased,  the  poles  of  the  system  move  to  the 
right  and  cross  over  into  the  right  half  plane 
indicating  that  the  system  is  no  longer  stable.  A 
Matlab  subroutine  to  generate  Figure  5  is 
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relatively  straightforward  to  write,  and  students 
should  be  able  to  easily  reproduce  the  figure. 

In  the  late  1 980’s,  Anderson  and  Spong  [3] 
attributed  the  instability  associated  with  time 
delay  to  a  non-passive  communication  block  and 
found  that  it  is  possible  to  stabilize  a  force 
reflecting  teleoperation  system  that  has  a  time 
delay  by  exploiting  scattering  theory.  Later, 
Niemeyer  and  Slotine  [10]  presented  the  wave 
variable  method  in  a  more  intuitive,  physically 
motivated  formalism  based  on  passivity  [10].  In 
the  wave  variable  method,  wave  variables  are 
used  in  place  of  more  conventional  power 
variables  such  as  velocity  and  force.  It  can  be 
shown  that  that  by  transforming  forces  and 
velocities  into  wave  variables  for  transmission 
along  the  communication  channel,  systems  could 
remain  stable  for  any  fixed  time  delay.  Several 
researchers  including  Munir  and  Book  [8]  have 
addressed  the  issue  of  stabilizing  systems  with 
variable  time  delay  by  incorporating  predictors 
into  the  system  along  with  a  conventional  wave 
variable  controller.  This  technique  was 
successfully  implemented  on  a  teleoperation 
system  that  used  the  Internet  (which  is  subject  to 
variable  time  delay)  as  the  communication 
channel. 

The  wave  variable  method  is  based  on 
passivity,  which  is  related  to  the  intuitive 
physical  concepts  of  power  and  energy.  The 
concept  of  passivity  provides  a  simple  and 
powerful  tool  for  analyzing  the  stability  of  a 
nonlinear  system.  It  is  based  on  the  intuitive 
notion  that  a  system  is  passive  if  it  absorbs  more 
energy  than  it  produces.  Space  limitations  do  not 
allow  us  to  provide  a  detailed  discussion  on 
passivity  theory  here.  Instead,  we  will  give  only 
a  quick  description  of  the  most  important  points 
of  passivity  and  scattering  theory  as  they  relate  to 
wave  variables  and  refer  the  interested  reader  to 
reference  [15],  which  provides  an  excellent 
introduction  to  a  wide  variety  of  nonlinear 
system  concepts  including  passivity  theory,  and 


reference  [9],  which  provides  more  details  on  the 
role  of  passivity  in  the  wave  variable  method. 

In  the  passivity  formulation,  one  defines  a 
power  quantity  Pin  entering  the  system  as  the 
scalar  product  between  the  input  vector  x  and 
output  vector  y  of  the  system.  Not  unlike  the 
relationship  between  Lyapunov  functions  and 
energy,  this  quantity  may  or  may  not  correspond 
to  a  physical  power.  Using  the  power  analogy, 
the  input  power  should  be  stored  or  dissipated  in 
the  system.  More  specifically,  let  Eslore  be  a 
lower  bounded  energy  storage  function  and  let 
Pdin  be  a  non-negative  power  dissipation 
function.  The  system  is  passive  if  it  satisfies 

Pln=jtE,tm  +  Pdl„.  0) 

If  the  power  dissipation  is  zero  for  all  time 
regardless  of  input  or  output,  the  system  is  said 
to  be  lossless;  however,  if  the  non-negative 
power  dissipation,  Pdiss ,  is  positive  during  any 
power  transfer  then  the  system  is  said  to  be 
dissipative.  Passivity  is  a  sufficient  condition  for 
stability. 

One  particularly  nice  feature  of  passivity  is  the 
important  property  that  if  multiple  passive 
elements  are  connected  together,  the  resulting 
system  is  also  passive.  This  ensures  that  a 
connection  of  passive  systems  remains  stable. 

The  stability  of  a  teleoperation  system  can  be 
proven  using  scattering  theory  [3].  Given  a 
transfer  function  g(j),  the  stability  of  the  system 
is  related  to  the  norm  of  the  scattering  matrix 

S(ja>)  =  [GUco)  -  I\[G(ja>)  +  /]"' .  (5) 

More  specifically,  for  the  system  to  be  passive, 
Ill'll  <1  where 

M  =  sup||S(./ctf)||2 .  (6) 

o 

The  norm  on  the  right  of  equation  (6)  is  the 
induced  2-norm,  i.e.,  \\s\\2=  ^m„{s'(j<a)S(j<»)) 
where  S\ja>)  denotes  the  complex  conjugate 
transpose  of  S(jco)  and  Xmn{A)  denotes  the 
largest  eigenvalue  of  the  Hermitian  matrix  A. 
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To  the  see  the  difficulty  caused  by  time  delay, 
suppose  that  the  communication  lines  in  Figure  1 
both  have  a  Tsec  time  delay.  The  transfer 
function  relationship  across  the  communication 
lines  is  then  given  by 


where 


Tm 


G{s)  = 


It  is  straightforward  to  show  that 
corresponding  scattering  operator  (5)  is 
_  .  \~  j  ian^coT)  sec (coT) 

-  sec (coT)  -  j  tan (coT) 


(7) 

(8) 

the 

(9) 


which  has  an  unbounded  norm  implying  that  the 
system  is  not  passive.  In  fact,  one  can  show  that 

[tan2(«7’)  +  sec2(ft)7')  }2sec(a)T)tar\(a>T) 

■  j2sec (euT)  tan(coT)  iar\l  (o)T) + sec2(w'/')J 


S'(jo})S(ja>)  = 


(10) 

■J 

which  has  eigenvalues  [tan(a>7)±sec(G>7)]  so 
that 


||S||  =  sup  |tan(n/T)|  +  |sec((u7,)|  =  oo .  (11) 

(o 

Once  again,  this  proves  that  the  scattering 
operator  is  unbounded  and  hence  the  system  is 
not  passive.  In  practice,  the  signals  entering  the 
communication  block  are  band-limited. 
However,  since  |tan(cu7)|+|sec(cu7)|  >  1  for  any 
co,  T>0  with  coT  i-  nit,  the  communication  block 
is  never  passive  for  any  range  of  frequencies. 
Furthermore,  we  would  expect  the  instability 
problems  to  be  most  severe  when  ooT-+it/2. 
Hence,  the  system  will  remain  stable  only  if  it 
exhibits  extremely  band-limited  behavior  [3]. 


Fortunately,  the  wave  variable  method  can  be 
implemented  to  stabilize  a  system.  This  is 
accomplished  by  using  wave  transformations  on 
both  the  master  and  slave  sides.  The  wave 
transformation  relations  for  the  single  degree-of- 
freedom  case  are  given  by 

u  (t)  =  u  At  — T) 

(12) 


v„(/)  =  vt(/-n. 


The  wave  transformations  for  the  left  wave 
junction  are  given  by 


Figure  6:  Wave  transformations  for  a  single 
degree-of-freedom  teleoperation  system. 


«.(0  = 


Vm(0  = 


•J2b 

■Jib 


(13) 


and  that  for  the  right  wave  junction  are  given  by 
bGAt)  +  Tc{t) 

Jib  (14) 

bGsAO-tJO 


«,(*)  = : 


v,(/)  =  - 


Jlb 


as  shown  in  Figure  6.  Although  the  strictly 
positive  parameter  b  can  be  chosen  arbitrarily,  it 
defines  a  characteristic  impedance  associated 
with  the  wave  variables  and  directly  affects  the 
system  behavior  [10],  The  transfer  system  across 
the  communication  line  now  becomes 
T_ 


Se 


\d  _ 


where 


Gw(*> 


b(l-e 


2e 


=  GAs) 

rJ. 

el!T) 

-le'T 

f2  *T 
-sT 

1  +  e~lsT 
l-e-2‘T 

(15) 


l+e~2U 

6tanh(jr) 

sechlsr) 


b(\  +  e~2sT) 
-  sech(j7) 
|tanh(jT) 


(16) 


This  can  be  easily  shown  by  applying  Mason’s 
Gain  Formula  to  the  block  diagram  shown  in 
Figure  6.  Since  sinh(/<9)  =  jsind,  cosh (jff)  =  cosd, 
and  tanh (jff)  -jiatcid,  it  follows  that 


GwUa>)  = 


(17) 


|  jb  tan  (coT)  -  sec(  coT) 
sec  {coT)  j  j  tan(«7’)j 

One  can  then  verily  directly  that  the  norm  of  the 
scattering  operator  S(jm)  =  [Gw{jo})-  i][Gw(jco)  +  i]-' 

is  ||S||  - 1 ,  which  means  that  the  system  is  passive 
and  lossless.  This  is  a  somewhat  tedious 
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Figure  7:  SIMULINK  block  diagram  of  a  bilateral  teleoperation  system  with  time  delay  and  the  wave 
variable  method. 


calculation;  however,  one  can  instead  note  that 
Gw(j<£>)  is  skew-Hermitian  and  then  use  the  fact 
that  the  scattering  operator  of  any  skew- 
Hermitian  matrix  is  unitary  and  hence  has  a  unit 
2-norm,  i.e.,  ||£(_/fi>)||  =  1  >  regardless  of  the  value 

of  co. 

While  the  previous  calculations  are  at  times 
tedious,  it  is  instructive  for  students  to  verify  the 
above  results.  Such  an  exercise  would  also  allow 
an  instructor  the  opportunity  to  reiterate  certain 
important  matrix  theory  concepts  such  as  matrix 
norms  and  unitary  matrices. 

To  incorporate  the  wave  variable  equations  into 
the  SIMULINK  model,  (13)  and  (14)  are 
rearranged  so  they  can  be  implemented  using 
simple  gain  blocks.  Figure  6  shows  this  simple 
configuration.  Once  the  wave  variable  method  is 
arranged  in  this  form,  it  is  easily  implemented  in 
SIMULINK  as  shown  in  Figure  7. 

Because  wave  variables  are  transmitted  across 
the  communication  lines,  the  system  is 
guaranteed  to  remain  stable  for  any  fixed  time 


delay.  As  can  be  seen  in  Figure  8,  the  system 
indeed  remains  stable,  and  the  force,  position, 
and  velocity  of  the  slave  clearly  follow  those  of 
the  master.  If  these  results  are  compared  to  the 
original  system  without  any  time  delay  present,  it 
can  be  seen  that  the  general  trends  are  preserved; 
however,  there  is  a  slight  degradation  in 
performance,  both  with  respect  to  peak  overshoot 
and  settling  time.  This  compromise  of  system 
performance  for  guaranteed  stability  is  a  natural 
part  of  control  system  design.  Overall,  even  with 
these  performance  issues  it  can  be  seen  that  the 
system  with  the  wave  variable  method  yields 
more  desirable  results  for  a  bilateral 
teleoperation  system  that  may  experience  time 
delay. 

A  nature  step  after  addressing  issues  associated 
with  constant  time  delay  is  to  consider 
teleoperation  systems  experiencing  varying  time 
delay.  This  occurs  for  example  when  the  Internet 
is  used  for  the  communication  link.  This  requires 
more  advanced  techniques  since  the  wave 
variable  method  by  itself  only  guarantees 
stability  in  the  presence  of  a  fixed  time  delay. 
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Figure  8:  Torque,  position,  and  velocity  data  for 
the  simulated  model  of  a  bilateral  teleoperation 
system  with  a  time  delay  given  by  2T=500msec, 
now  with  the  wave  variable  method 
implemented. 

Schemes  using  additional  elements  such  as 
observers  and  Kalman  filters  as  estimators  have 
successfully  stabilized  teleoperation  system 
using  the  Internet  [8];  however,  the  issue  of 
system  performance  is  still  an  ongoing  research 
problem. 

IV.  Conclusions 

In  summary,  we  have  described  the  Matlab  and 
SIMULINK  tools  necessary  to  simulate  and 
control  a  bilateral  teleoperation  system.  Matlab 
was  used  to  set  up  the  parameters  of  the  system, 
and  SIMULINK  was  used  to  set  up  the 
connections  between  the  master  and  slave 
models  and  to  simulate  the  overall  system.  After 
simulating  a  simple  teleoperation  system 
operating  under  ideal  conditions,  time  delay  was 
introduced  into  the  communication  link  of  the 
model.  It  was  shown  that  even  a  modest  amount 
of  time  delay  can  be  detrimental  to  a 
teleoperation  system  using  force  feedback.  A 
relatively  new  technique  for  handling  time  delay 
called  the  wave  variable  method  was  then 
implemented  to  stabilize  the  teleoperation  system 
with  time  delay.  Other  possible  problems  to 
include  in  a  classroom  project  are  noise 
compensation  and  variable  time  delay.  Although 


this  article  focused  on  single  degree-of-freedom 
systems,  the  system  described  here  can  be 
extended  to  a  multiple  degree-of-freedom  system 
[1,2]. 

The  material  presented  in  this  article  should  be 
appropriate  for  inclusion  in  an  undergraduate  or 
first  year  graduate  course  in  control  systems.  The 
material  only  requires  access  to  a  PC  with 
computer  software  such  as  Matlab  with 
SIMULINK. 
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